Optical and magneto-optical far-infrared properties of bilayer graphene 
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We analyze the spectroscopic features of bilayer graphene determined by the formation of pairs 
of low-energy and split bands in this material. We show that the inter-Landau-level absorption 
spectrum in bilayer graphene at high magnetic field is much denser in the far-infrared range than 
that in monolayer material, and that the polarization dependence of its lowest-energy peak can be 
used to test the form of the bilayer ground state in the quantum Hall-effect regime. 



Monolayer and bilayer graphenes^ 2 -! 3 - are gapless two- 
dimensional (2D) semiconductors^ 5 -. When used in 
transistor-type devices, the density and type (n or p) 
of carriers in them can be controlled using the underly- 
ing gate 1 , which has been exploited in the recent stud- 
ies of the quantum Hall effect (QHE) in such struc- 
tures. The sequencing of QHE plateaus in graphenei* 2 -' 3 - 
revealed the peculiar properties of the charge carriers 
in this material: Dirac-type chiral electrons with Berry 
phase 7r in monolayers^, and the Berry phase 27r chi- 
ral quasiparticles 5 with doubled degeneracy of the zero- 
energy Landau level (LL) in bilayers. 

The electromagnetic (EM) field absorption in graphene 
at zero magnetic field has already been studied££&i£. 
While the dc conductivity of monolayer graphene in- 
creases linearly with the carrier densitjiii, the real part 
of its high-frequency conduct ivitji^ is independent of 
the electron density in a wide spectral range above the 
threshold hio > 2|ep|, which determines a featureless ab- 
sorption coefficient g\ — ire 2 /he. In contrast, the bilayer 
absorption coefficient, 



fj2 



(2ire 2 /Hc)f 2 (u), 



(1) 



reflects the presence and dispersion of two pairs of bands 
in this material 5 ^ 12 ' 13 . In this paper, we analyse far- 
infrared (FIR) magneto-optical properties due to the 
inter-LL^ i 14 ! 15 transitions in bilayer graphene and com- 
pare them to those in the monolayer material. The low- 
energy part of the absorption spectrum in a bilayer sub- 
jected to a strong perpendicular magnetic field is much 
denser than that of a monolayer, reflecting the parabolic 
dispersion of its low-energy bands, with the lowest FIR 
absorption peak formed by transitions involving one of 
two degenerate LLs at zero energy. Although the high- 
energy part of the bilayer absorption spectrum is not sen- 
sitive to the polarisation of FIR irradiation, below we 
show that its lowest peak may be strongly polarised re- 
flecting the occupancy of the degerate LLs in a bilayer 
determined by the form of its ground state in the QHE 
regime. 

The electronic Fermi line in graphene surrounds the 
corners^ K± of the hexagonal Brillouin zone^ (where 
we set e = 0). In monolayer graphene, quasiparticles 
near the centres of valleys K± can be described by 4- 
component Bloch functions ip = [</>k +j 4, <Pk. + b, ^k_b, 
</>k_a] and the Hamiltonian Hi — vll z ® <x ■ p, where 
c = ipxiVy) are Pauli matrices acting in the space of 



electronic amplitudes on the two crystalline sublattices 
(A and B) and H z is the diagonal Pauli matrix in the 
valley space. Momentum p = — iKV — -A is calculated 
with respect to the center of the corresponding valley and 
VxA = Bl 2 . 

Bilayer graphene is composed of two coupled mono- 
layers (with sublattices A, B and A, B in the bottom 
and top layers respectively) arranged according to Bcrnal 
stacking^: sites B of the honeycomb lattice in the bottom 
layer lie below A of the top layer. It also has a hexago- 
nal Brillouin zone with two inequivalent valleys, but car- 
ries twice the number of electronic dispersion branches. 
The latter can be found using the nearest-neighbour 
hopping Hamiltonian 5 acting in the space of sublattice 

States </> K+ B> ^K+A' <^K + s; </> K _.B> </>K_B, 



Hi 



a t p vll z <g! 

<t p 7i n 



<X p 



(2) 



Here, v is determined by the AB (AB) intra-layer hop- 
ping, 71 is the strongest inter-layer AB hopping clement, 
and 7J3 <C v is due to a weak direct AB hop, and 't' stands 
for transposition. Equation ^ determines two types of 
branches in the electronic spectrum of graphene 5 -^ ob- 
served in the recent ARPES studies^ 3 : Split-bands, 



ef = ±( 7 i/2)U/l + 4 U V/7i 2 + l], 



(3) 



formed by symmetric and antisymmetric states based 
upon AB sublattices; and two gapless branches, ef. For 
e > el s 371 (f 3 AO 2 , the dispersion in gapless branches 
can be approximated by 



±(71/2)^1 + Av 2 p 2 /j 2 -1]. 



(4) 



For |e| < jji, gapless bands are formed by states from 
sublattices A and B described using the 4-component 
Bloch functions \ = [<f>K + A, <fi K+ B> ^k_S' 0k_a] and 
the reduced low-energy Hamiltonian 5 -^, 



H 2 = - 
Sh w = v 3 U 



1 - f3p 2 
(<r • p)cr x {cr ■ p) + Sh w H ; 

2m2 m 2 



er* • p; m 2 = 7i/2f 2 



< 1. 



(5) 



Here m 2 « 0.05m e and the last term in H 2 takes into 
account the weak AA and BB intra-layer hopping 18 . 
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In a 2D electron gas with conductivity <j(lu) much less 
than c/2-k, absorption of an EM field E w = £Ee~ tuJt 
with polarisation t (£® = [l x — il y ]/y/2 for right- and 
£q = [l x + i\ y ]/y/2 for left-hand circularly polarised 
light, respectively arriving along the direction antiparallel 
to a magnetic field) can be characterised by the absorp- 
tion coefficient g = EiE*<Jij(uj) / S: the ratio between 
Joule heating and the energy flux S = cE x H/47T = 
— S\ z transported by the EM field. Using the Keldysh 
technique, we express 



8e< 



5i 



Fde, 



-Tr 



where v = d p H is the velocity operator, Tr includes the 
summation both over the sublattice indices "tr" and over 
single-particle orbital states, N is the normalisation area 
of the sample and F = n-p(e) - np(e + u>) takes into 
account the occupancy of the initial and final states. Here 
we have included spin and valley degeneracy. 

For a 2D gas in a zero magnetic field, the elec- 
tron states are weakly scattered plane waves. Us- 
ing the plane wave basis and the matrix form of the 
high energy bilayer Hamiltonian TC 2 , we express the re- 
tarded/advanced Greens functions of electrons in the 



bilayer as G R / A (p, e) 



e± \ihr- 1 



H 2 (P) 



and 



Tr = j d 2 p / 2 ^La tr, neglecting the renormalisation of the 
current operator by vertex corrections at lot 3> 1 and the 
momentum transfer from light (since v/c ~ 3 x 10 -3 ). 
This reproduces the constant absorption coefficient g\ = 
Tte 2 /hc so that f\ — | 7 ' 9 m monolayer graphene and 
yields the following expression for the absorption coeffi- 
cient of bilayer graphene for light polarised in the plane 
of its sheeti^: 



92 



he 



71 71 

0(fi-l) (ft -2)0(0-2) 



2(0 + 1) 



Q 2 



2(n - 1) 



(6) 
(7) 



Here 6(x < 0) = and 9{x > 0) = 1. This result agrees 
with the calculation by J. Nilsson et al£- taken in the 
clean limit and at T = 0. The frequency dependence 20 
of the bilayer optical absorption is illustrated in FigJU 
where an additional structure in the vicinity of Hlj — 
71 (71 w 0.4eV (Refs^ ancU 3 -) is due to the electron- 
hole excitation between the low-energy band ef- and the 
split band e~. For the higher photon energies, frui > 
271, the frequency dependence almost saturates at / « 
1. Over the entire spectral interval shown in FigQ] the 
absorption coefficient for the left- and right-handed light 
are the same, so that Eq. (JJJ is applicable 19 to light 
linearly polarised in the graphene plane. 

High-field FIR magneto-optics of graphene. In a mag- 
netic field, the continuous zero-field spectrum, [Eq. {7}] 
splits into lines determined by transitions between Lan- 
dau levels. The LLs can be studied after rewriting Hi 
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FIG. 1: Absorption coefficient of bilayer and monolayer 
graphene in the optical range of frequencies. Insets illustrate 
the quasiparticle dispersion branches in the vicinity of cf and 
possible optical transitions. 



and H 2 in terms of descending, n = p x + ip y and raising, 
71^ = p x — ip y operators in the basis of Landau functions 
tp n >o, leading to 



V 



7T+ 
7T 



The resulting monolayer spectrum^ contains 4-fold de- 
generate (2x spin and 2x valley index) states: one 
at e = with tp 0K+ = [ip , 0,0,0] and ip 0K _ = 
[0,0,ipo,0], and pairs of levels e n ± — ±hv\^ l \/2n 
hc/eB is the magnetic length) with ip n ±,K + = 
-1,0,0] and ip n ± l K- = ^[0, 0, ip n) =R¥>n-l]- 



-±=[(p n) ±icp n 

We neglect the electron spin splitting and use the index 
a = ± in e na and ipna.K for the conduction (+) and 
valence (— ) band LLs. 

The bilayer spectrum features^ a group of eight states 
with Eq ~ £1 ~ with wave functions xo.k + = 
0,0,0], xo,K- = [0,0,</J ,0] and X i,k+ = [(pi, 0,0,0], 
Xi,K_ = [0)0, <fi,0] (4x2, due to spin degeneracy), and 
a ladder of al most equ idistant 4-fold degenerate levels 
e n ± = ±huj c yjn{n — I) (2x spin and 2x the valley in- 
dex) with hu>c = h 2 jm 2 \ 2 B and wave functions Xn±,K + = 
^bn,±¥>n-2,0,0], Xn±,K- = ^[0,0,</J n ,±^„_ 2 ]. The 
latter result can be found from the analysis of the first 
term in the low-energy-band Hamiltonian H 2 , which 
dominates for |e| < ^71 and high magnetic fields such 
that Xg 1 > 71U3/U 2 . Numerical diagonalisation of the 
full H 2 and 7i 2 shows^ that this degeneracy is not lifted 
by the warping term 5H W , and also the above-described 
grouping of LLs in bilayer graphene was confirmed in 
the recent QHE measurements 3 -. It has been noticed^ 
that degenerate levels with n = {xo,k+ 1 Xo,K- ) an d 
n = 1 (xi,k+, Xi,k ) can be weakly split by the second- 
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neighbour hopping term (3p /m?, to Eq — ^ftiu c and 
Ei = ^-hjj c , where (3 -C 1— • However, it is more likely 
that, at a high magnetic field, electrons in an ideally clean 
bilayer with filling factor —4 < v < 4 would form a corre- 
lated ground state in which the occupancy of n — and 
n = 1 LLs would be determined by the electron-electron 
interaction. The correlated ground state may be particu- 
larly interesting in a bilayer with v = 0. One can envisage 
2D electrons forming a ferromagnetic QHE state in which 
both n = and n = 1 states are half-filled (like in a con- 
vensional QHE liquid at the filling factors corresponding 
to a half-filled spin-degenerate LL). Alternatively, there 
can be an antiferromagnetic state with one fully occu- 
pied and one empty LL. Below we show that these two 
ground states can be distinguished experimentally on the 
basis of magneto-absorption spectra measured in circu- 
larly polarised FIR light. 

Figure [Ha) illustrates the selection rules for the single- 
particle inter-LL transitions in bilayer and monolayer 
graphene. In both cases, photons with in-plane^ po- 
larisation £q are absorbed via transitions where the LL 
index changes from n to n — 1 (n > 1), whereas ab- 
sorption of £® photons happens via transitions from the 
LL n to n + 1 (n > 0) . Assuming the same broadening 
hr^ 1 of all LLs— , we arrive at the following magneto- 
absorption spectra for monolayer^ and bilayer graphene 
(for e L < hu < ±7!): 



ffi/e/ D \ 27re 2 j.e/G/ D \ 



(8) 



2h 



ff 



E 



^ (te _ a ^T^l + a'VH) 2 + 1 



2b„ / _ x 



Kit 



V 

n> 
act 



o^(£-aV»+T+«V») +1 



-fn-l,a)(n-l) 



f 2 e = E 



a\/ n 2 — n— a' yj (n— l)(n— 2) 



2 - aVn 2 - n+a'VTn - - 2) 

/ 2 e = E 



7tcj c t Qi/n 2 +^— a' y/n 2 — r 



n>l 
a a 



if—- a\/n 2 +7i + a'Vn 2 - n) + 



Here wi = V^hvXg 1 and Sw c = h 2 lm?.\\ are the char- 
acteristic energy scales for the LL spectra in monolayer 
(J = l;e„ Q ) and bilayer (J = 2;e na ) graphene, respec- 
tively, and a, a' = ± determine whether the correspond- 
ing state belongs to the conduction (+) or valence (— ) 
band. Also, v n a are the filling factors of the correspond- 
ing LLs, and bo — 1, b n >i = l/\/2 for a monolayer and 
co,i = 1, c„>2 = l/v2 for the bilayer. 

One can test the selection rules shown in Fig. [2^a) and 
the actual polarisation of these transitions using gated 
graphene structures^ 2 -!-?. By filling the monolayer sheet 
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FIG. 2: (a) Allowed inter-LL transitions without trigonal 
warping effects. Dashed and dash-dot lines indicate transi- 
tions in £® and Iq polarisations respectively, (b) Monolayer 
(top) and bilayer (bottom) FIR absorption spectra in £9 and 
Iq polarisations for B = 10T and filling factor v = 0. Dashed 
and solid lines describe absorption by ferro- and antiferro- 
magnetic states of the v = bilayer. (c) Weak field, B = IT 
magnetoabsorption in bilayer graphene with v = calculated 
with ( 1)3 = 0.2u) and without (113 = 0, shaded) warping term 
in the Hamiltonian. 



with electrons up to v — 2 (a completely filled n — 
LL), one would suppress the intensity of the lowest ab- 
sorption peak in £q polarisation and increase the size of 
the £® peak. By depleting the monolayer to v = —2 state 
(emptying the n = LL) one would achieve the opposite 
effect. Similarly, in a bilayer with a completely filled pair 
of n = 0, 1 LLs which takes place at v = 4, light with 
can be absorbed only in the £® polarisation. 
By depleting the bilayer down to v — —4 one could sup- 
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press the absorption in this line in €q polarisation while 
retaining £q absorption. 

The denser magneto-absorption spectrum in bilayer 
graphene (as compared to the monolayer) is one of its 
features illustrated in Fig|2] (b) . Another feature is that 
the intensity of the lowest absorption peak measured in 
the £® or £q polarisations at the filling factor v = 
may differ for different ground states of the bilayer. To 
illustrate this possibility, let us compare the absorption 
spectra for two model ground states of the bilayer. 

The first can be attributed to a bilayer with such a 
large single-particle splitting between the n = and 
n = 1 states that one of these levels is full and the 
other empty, even in the presence of the electron-electron 
interaction. This may be due either to the inter-layer 
assymmetry£ii£ or caused by the AA/BB intralayer hop- 
ping in the last term of H% , Equation {1} . Below we refer 
to such a QHE state as being antiferromagnetic, stressing 
that this state is not spin-polarised. 

An alternative form of the ground state can be at- 
tributed to the case of a negligible splitting between the 
n = and n = 1 LLs. In the latter situation the electron- 
electron repulsion may lead to the ferromagnetic align- 
ment of electron spins and the formation a ferromagnetic 
QHE state with one spin component of each LL being 
completely full and the other completely empty. 

Since transitions from/to the n — LL with £o w 
to/from the states e n ± with n > 2 are forbidden, the 
intensity and polarisation of the peaks at u = \Z2ui c are 
determined by transitions from/to n = 1 LL (also with 
Ex w 0) and directly reflect the occupancy of this state. 
If the bilayer ground state is ferromagnetic, with a half- 
filled n = 1 LL, the absorption peak at ui = V2uj c will 
have the same intensity in both I® and £q polarisations. 

In contrast, absorption by a v = bilayer with an- 
tiferromagnetic ground state would contain the line at 
u> = \/2w c only in one polarisation: in £q if n = 1 LL 
is empty (fully occupied n — LL) and in £^ if n = 1 
LL is full. For comparison, the lowest peak in the spec- 
trum of a monolayer with v = appears in both polar- 
isations, since both transitions to (ei_ — > eo) and from 
(eo — > ei+) half- filled n = monolayer LL are possible 
All higher-energy absorption peaks which involve tran- 
sitions between filled valence band states (a = — ) and 
empty states in the conduction band (a = +) have equal 
strengths in both polarisations, which reflect effectively 
the inter-band nature of these transitions. Also, note 
that a weak transition between n = and 1 LLs with a 
low frequency (3lo c (microwave range) is possible, to the 
measure of a small (3 and depending on the LL filling. 

We briefly turn our attention to the effect of the 
electron-electron interaction on the absorption. Since we 
have specified that we study the spectrum for exactly 
filled levels we can construct magneto-exciton operators 
for each case^ 2 -. These operators correspond to collective 
excitations of electrons from fully filled to completely 
empty LLs. We calculated the shift in the dispersion 
of such magneoexcitons caused by these interactions to 



the first order in the interaction constant e 2 /(e r Asfio; c ) 
- where e r is the dielectric constant of graphene - and 
find that the Coulomb interaction has no effect on these 
excitations when their in-plane wave vector k is zero. 
This result agrees with recent studies of the mangetoex- 
citon dispersion in the integer QHE states in monolayer 
graphene 2 ^ and our previous work in the context of semi- 
conductor heterostructures probed with surface acoustic 
waves 2 ^. Therefore in these cases we can disregard the 
electron-electron interactions and say that our conclu- 
sions concerning the polarisation properties of FIR ab- 
sorption apply to all magnetoexcitons which are sym- 
metric in the spin and valley indices. 

The analysis mentioned above was performed in the 
lowest order of pertubation theory in the interaction. Its 
result does not imply a full Kohn's theorem as known for 
electron-electron interactions in a simple parabolic band. 
In graphene, the chirality of the carriers is similar to 
the spin-orbit coupling in III-V semiconductors in that it 
plays the role of a non-parabolicity in the band structure 
which is known to violate Kohn's theorem in conventional 
semiconductors. This means that the conclusions for the 
lowest order corrections due to the interaction may not 
stand for higher orders of pertubation theory. This issue 
will be the subject of a seperate investigation 25 . 

Finally, the warping term Sh w in Eq. {5| mixes each 
state Xna.K with states X(n±3)a,K- This generates weak 
transitions (with the coefficient <5<?2 ~ ^31x12X3 / ti) 2 32) 
between states which are 2 and 4 levels apart. Although 
such transitions are negligibly weak at high fields where 
the first term in H2 is dominant, they become relevant at 
weak fields, where A^ 1 < liv^/v 2 , and the low-frequency 
absorption spectrum of a bilayer aquires an additional 
structure. In Figf5)(c) we compare the absorption spectra 
m a v = bilayer. 2 ! at B = IT calculated numerically 
for U3 = and V3 = 0.2v. 

In conclusion, we have described peculiar optical and 
FIR magneto-optical properties of bilayer graphene and 
compared them with those of the monolayer material. 
The zero field optical absorption spectrum shown in Fig. 
Q] reflects the existance of four bands in the electronic 
spectrum of this material whereas the FIR magneto- 
optical spectrum of bilayer graphene, which is much 
denser than that of a monolayer, reflects the parabolic 
dispersion of its low-energy bands e f . Due to the inter- 
band character of the inter-LL transitions in graphene 
with low carrier density, the magneto-optical spectrum 
in a monolayer and the higher energy part of the bilayer 
spectrum (w > 3oj c ) do not depend on the FIR light po- 
larisation. However, the absorption peak at lu — ^/2uj c 
may appear differently in and £ e polarised light, de- 
pending on the occupancy of the two zero-energy LL 
(n = and 1) in bilayer graphene^i 5 -. In such spectro- 
scopic studies, the weakness of (72 can be overcome once 
one uses the 2D electrons in graphene in the QHE regime 
for the transport detection of FIR absorption. 
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